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Abstract 

We use numerical models, supported by our laboratory data, to predict the 
dust densities of ejecta outflux at any altitude within the Hill spheres of 
Europa and Ganymede. The ejecta are created by micrometeoroid bombard- 
ment and five different dust populations are investigated as sources of dust 
around the moons. The impacting dust flux (influx) causes the ejection of 
a certain amount of surface material (outflux). The outflux populates the 
space around the moons, where a part of the ejecta escapes and the rest falls 
back to the surface. These models were validated against existing Galileo 
DDS (Dust Detector System) data collected during Europa and Ganymede 
flybys. Uncertainties of the input parameters and their effects on the model 
outcome are also included. The results of this model are important for fu- 
ture missions to Europa and Ganymede, such as JUICE (JUpiter ICy moon 
Explorer), recently selected as ESA's next large space mission to be launched 
in 2022. 
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1. Introduction. 



Micrometeoroids (solid micron-sized dust particles) are a common con- 
stituent of the Solar System. They can easily reach the surfaces of atmo- 
sphereless bodies and, upon impact, cause the ejection of surface material 
into the surrounding space. The ejected dust fragments populate the space 
around the host bodies, where a part of the ejecta escapes and the rest even- 
tually falls back to the surface. This process can be numerically modelled. 
In this paper we investigate dust around Europa and Ganymede created by 
micrometeoroid bombardment. 
Our study is built upon the models developed by 



Kriiger et al. 



Krivov et al. 



(120031 ) supported by our impact experiments (jMiljkovic et al. 



(2003) 



ana 



201 ll ) and hydrocode impact modelling to reduce the number of variables 



in the model. The impact experiments were made using the light gas gun 



(LGG) at the Open 



( jMiljkovic et al. 



2011 



Jniversity's 



Patel et al. 



Hypervelocity Impact 



2010 



McDonnell 



2006 



(HVI) laboratory 



Taylor et al 



20061) 



A series of high velocity impacts were made at 2 km s 1 using 1 mm diam- 
eter stainless steel balls as projectiles in to pure water ice and sulphate hy- 



drated mineral targets. The ejecta size ( jMiljkovic et al. 



20111 ) and velocity 



distributions were measured and subsequently modelled us ing ANSYS AU 



TODYN finite eleme nt hydro-dynamic shock physics code (jMiljkovid . 



Pierazzo et al. 



2010). 



2010 : 



The dust cloud model presented in this paper characterizes the dust envi- 
ronment (size, density, flux and velocity distribution of such dust) around 



Europa or Ganymede, predicts the dust densities of the ejecta outfhrx at any 
altitude within Europa's and Ganymede's Hill spheres (radius of 13300 km 
or 8.5 Re, where Re = 1565 km for Europa and 32000 km or 12 Rq, Rq = 
2631 km for Ganymede) and can evaluate the dust density at any altitude as 
a function of the size distribution of the dust. We choose a discrete selection 
of altitudes and dust masses to give a representative set of results. 
Our results are important for future space missions to Jupiter System that 
carry a dust detector onboard. This study can be further applied to esti- 
mate the dust counts into a dust detector in orbit around Europa and/or 
Ganymede. A dust detector has been proposed as part of a payload for a 
space mission to Europa and Ganymede. Initially named Laplace in 2007, 
it was renamed the Euro pa- Jupiter System Mission after ESA and NASA 



joined proposals in 2009 ( 



Blanc et al. 



20091 ) for a major mission to Jupiter 



System. In 2011, EJSM was reformed again into an ESA- led mission to 
Ganymede with flybys to Europa and C alisto, named JUICE (JUpiter ICy 



moon Explorer) ({Dougherty et al 



201 ll ). JUICE has recently been selected 



as ESA's next large mission for launch in 2022. If a dust detector is included 
in the payload, an in-situ analysis of the dust that surrounds Europa and 
Ganymede will be possible, which could provide information abo ut the sur- 



face, as its composition should be "written" in the detected dust (iMiljkovid . 



201lh . 



The proposed dust detector should not only be capable of determining the 
density of dust in the cloud, but may provide chemical analysis of cap- 
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tured dust (jMiljkovic et all 120081 ). as was the case with the analysis of 
the Jovian Dust Strea m particles by Cassini's CDA (Cosmic Dust Analyser) 



(jPostberg et al. 



20061 ) . Chemical abundance maps of Europa and Ganymede 
even at low spatial mapping resolution show that water ice is non-uniformly 
distributed over the surfaces of Europa and Ganymede. The Galileo NIMS 
(Near Infra-red Mapping Spectrometer) spectra identified t hese impurities 



as hy drated minerals, sulphates and possibly hydrocarbons (IMcCord et al. 



19981 ). As the surface material is ejected by micrometeoroid bombardment, 
it can be expected that the dust particles around Europa will be composed 
of water ice, sulphur salts an d their decompo sition products, including any 



potential organic compounds ( jMiljkovic 



20111 ) . It should be emphasized that 



a dust detector has never visited Europa or Ganymede at such a close orbital 
distance or spent a longer time than in flybys. No chemical analysis of the 
dust in Jupiter System has been made yet, apart from Cassini's Cosmic Dust 
Analyser (CP A) measurements of th e Io stream dust at more than 1 AU away 



from Jupiter (jPostberg et al 



2006|) 



1.1. Comparison of the present dust cloud model with previous models. 

I 



In previous work, iKrivov et al. 



( 120031 ) developed a spherically symmet- 



ric case for an atmosphereless body with applications to Ganymede and 
the Saturnian satellites to help the interpretation of Ca ssini's CDA (Cos- 



mic Dust Analyser) measurements; iKriiger et al. 



(I200Q[ ) developed a dust 



cloud model around Ganymede to explain the Galileo DDS measurements; 
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Sremcevic et al. 



( 120031 ) investigat ed the asymmetry of th e dust clouds around 



Galilean and Saturnian satellites; 



Sremcevic et al. 



(120051 ) compared 



;he Galileo 



data w ith their models for Europa, Ganymede and Callisto, whereas 



Kriiger et al. 



(120031 ) investigated the dust clouds around all four Galilean sate llites. 



he d ifferences between the dust model presented here and that of lKriiger et al. 



(120031 ) are the following: their slope of the cumulative ejecta mass dis- 



tribution was approximated, whereas in our model our slope was derive d 



directly from the im pact experiments published in 



Kriiger et al. 



Miljkovic et al. 



(2011). 



(120031 ) reported that dust have mean mass of 10 _11 g around 
Europa and 10 _13 g around Ganymede, whereas in our model the ejecta frag- 
ment size distributions were calculat ed, and were in rang e between 10 _15 g- 



10 g. Ejecta speed distributions in 



Kriiger et al. 



( 120031 ) were represented 



as the distribution of ejecta material having speeds higher than a certain 
speed, that is dependent on the minimum ejecta speed of fragments and the 
power-law slope of the distribution, which were unknown variables fitted to 
match the Galileo data, whereas in our model, a size-velocity relation was 
applied to all the ejected fragments in order to have more precise outffux 
and spatial densities of ejected dust. We primarily focus on the short-lived, 
bound ejecta at a distance from the surface at which a spacecraft may orbit. 
Therefore, any asymmetry effects in the spatial density of ejected fragments 



caused by Europa's or 



Sremcevic et al. 



}ital motion can be excluded (such were considered by 



(120031 )) as well as the charging of dust fragments. 
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2. Micrometeoroid influx into Europa's and Ganymede's surface. 



In the Jovian system, there are five distinct sources of dust around Europa 
and Ganymede. These are: (i) asteroidal and (ii) halo dust populations, 

ation disting uished 



19931 ): (iii) 



as part of the interplanetary dust particle (IDP) popu 
by their location and not necessarily by their origin ((Divine, 
interstellar dust (ISD) that originate from beyond the Solar System; (iv) the 
Io stream and (v) ring dust that from the Jovian system itself. 

2.1. Influx into Europa's and Ganymede's surface from the Solar System and 
interstellar region. 



Divind (119931 ) created a phase density model to predict micrometeoroid 
fluxes at different distances from the Sun. According to this model, there 
are five distinct interplanetary dust populations, out of which only two (as- 
teroidal and halo) are present at Jupiter' s orbital (heliocentric) distance, r. 
The cumulative number of IDPs per unit volume (spatial density, Nm ) whose 
mass exceeds m can be calculated as a function of inclination, i, represented 
by elliptical latitu de A, eccentric ity, e, and the IDP mass represented by H m 



is shown in Eq. [T] (IDivind . 



1993|). 



N, 



Al 



71 



H m dm 



7T-|A| 



Pi sin i di 



n Je x V 6 ~ 6 X J\X\ ^(COSA) 2 - (cOSi) 2 



(1) 



The mass distribution of micrometeoroids H m is independent of the posi- 
tion and velocity of dust particles in the Solar System, the cumulative mass 
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distribution of dust particles Ni is a function of radial distance from the 



Sun; x = s i n_ ( r i/ r ) an d e x = (r — r\)/{r + ri), where r\ is the perihe- 
lion distance, r=5.2 AU and A = 0° are Jovian heliocentric distance and 
the equatorial latitude, p e and Pi are normalized eccentricity and l atitude 
distrib utions. Discrete values for Aq, H m , p e and pi were taken from 



Divine 



( 119931 ) for the respective asteroidal and halo populations and integrated in 
Eq. [TJ H m was integrated over mass bins (Am), in order to transform the 
continuous mass distribution into a binned one (Eq. [2]). 



"m+Am 

Nm = const I H m dm 



(2) 



In newer meteoroid codes, the product of functions Nixp e xpi is replaced by 
a single vector functi on, which proyides a more accurate and detailed dust 



flux calculation (e.g. 



Dikarev et al 



( 120051 )). However, due to lack of obser- 



vational data at 5 AU from t he Sun, and s ince later meteo r oid m odels were 
based on 



Divinel (I1993h fe. g. 



sider the model by 



Jehnl d2000D 



Staubach et al.l ( 119971 )). we con 



Divinel (119931 ) to be satisfactory for a preliminary dust 



flux calculation at Jupiter's distance for dust coming from asteroidal and halo 
sources. It should be noted that we have neglected Solar radiation pressure 
due to Jupiter's large distan ce from the Sun. 



Colwell and Horanyil ( 119961 ) calculated that at 100 Rj away from Jupiter, 
the Oort Cloud dust (highly inclined and eccentric dust associated with Di- 
vine's halo population) shown as triangles in Fig. [Hand planetary dust (low 
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inclination and low eccentricity orbits, associated with Divine's asteroidal 
population) shown in squares in Fig. (TJ move at a mean speed of 23.6 kms~ l 
and 6.6 fcms -1 , respectively. 

The interstellar dust (ISD) can be observed far above the equatorial plane. 
The Ulysses spacecraft monitored the ISD activity at high ecliptic latitudes 
between 3 and 5 AU from t he Sun, which was far away from contamination 



by IPD 



Griin et al. 



kms 1 ( jLandgraf et al. 



19971) . I SD penetrate s the solar system at about 26 



2000 



Kriiger et al. 



20071 ). Fig. H shows the ISD 



flux data taken by Ulysses spacecraft (that measured grain masses between 
10 -11 g and 10~ 7 g) and ground based radar meteor observations made by 
AMOR (Advanced Meteor Orbit Radar) facility, that measured the flux of 
ISD dust large r than 10~ 7 g. AMOR data in Fig. [TJ is shown as a tail with 



a slope of -1.1 (jLandgraf et al. 



20001 ) . Fig. [JJalso shows that the ISD flux is 
lower than the IDP flux. Both the IDP and ISD densities were approximated 
to 2.5 gcm~ 3 . 

All grains approaching Jupiter are exposed to the Jovian magnetospheric 
plasma and magnetic field with which they interact. The effect of the 
Lorentz force on grain motion depends on the grain charge to mass ratio 



( iHoranyi et al. 



1993af ). The smaller the grain, the stronger is the perturba- 
tion, for a given charge. In order to account for the effects of the Jovian mag- 
netosphere, the influx dust masses in the range 10 _14 -10~ 10 g were multiplied 
by a factor of 2. 65, which is an averag e d mag nification value from the calcu- 



lations made by 



Colwell and Horanyi 



(119961 ) . Jovian gravitational focusing 
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increases the speed of approaching micrometeoroids and therefore the spa- 
tial density of micrometeoroids becomes higher closer to the planet. Grains 
smaller than 10~ 15 g were excluded from our model as they are likely to be 



come deflected from entering Jupiter's magnetic field ( IColwell and Hor anvil . 
19961 ) . On grains larger than 10~ 10 g gravitational effects are more dominant. 



The gravitationally focused influxes at Europa' s and Ganymede's distances 



from Jupiter, F in , were calculated using Eq. [3] (ISpahn et al. 



20061 ). 



F in /F™ = 0.5VT+TIVT+T + 1 + W1 + T — (i? p /a) 2 (l + Ta/R p )} (3) 



F££ = Nmv°^ p is the non-focused micrometeoroid flux at Jupiter's heliocen- 



tric distance, for all different dust populations. T = 2GM p /(a(v, 



oo \2\ „.oo 
imp J )i u imp 



is the micrometeoroid non-focused velocity, Vi mp is the micrometeoroid grav- 
itationally focused velocity, G is the gravitational constant, M p is Jupiter's 
mass, a is the mean semi-major axis of Europa's or Ganymede's orbit and 
R p is Jupiter's radius. 

2.2. The influx onto Europa and Ganymede sourced within the Jovian system. 

There are also dust populations that originate from inside the Jovian sys- 
tem and intersect the orbits of Europa and Ganymede (Fig. [T]). Io stream 



dust (5-15 nm in size, assu ming the dust is spherical (IKriiger et al 



1.5-2.16 gem 3 in density ( IKriiger et al 



2004; 



Postberg et a 



spatia lly averaged dust density between 10 3 and 10 8 m 



(120041 )) and ring dust, 500-1000 km 3 in spatial dust density (IKriiger et al. 



200J), 



2006 ) and 



(Krii ger et al 
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2004 ; iKrivov et al.l . |2002| ). 0.6-2 /im in size, assuming sphe rical shape and 



2004 ) also contribute 



approximately 2 gcm~ 3 in bulk density ( iKriiger et al. . 
to the bombardment of Europa and Ganymede. Io dust quickly becomes 
charged and its mo tion is then influen ced by the Jovian and interplane- 



tary magnetic field (IFlandes et al 



201ll ). forming a so-called balle rina skirt 



shaped dust stream and moving at the co-rotational plasma velocity ( jHoranyi et al. 



1993bl ). Therefore it can be assumed that the Io stream dust approaches Eu- 
ropa at a speed of about 120 kins' 1 and Ganymede at about 190 kms' 1 . 
The dust rings are composed of dust grains moving in prograde and ret- 
rograde direction. Prograde dust is likely to have been ejected from the 
surfaces of the Galilean moons and assumed to be four times more abundant 



( iThiessenhusen et al. 



20001 ) than retrograde dust. Retrograde dust is likely 



to be populated by captured IDPs. It is not known if non-gravitationa l 



20021 ). 



processes affect the motion of the dust in the rings ( IKrivov et al. 
Europa and Ganymede move together with the prograde dust ring at ap- 
proximately 14 kms~ x and 11 kms -1 , respectively. Impacts onto Europa's 
and Ganymede's surfaces made by retrograde dust can be considered to oc- 
cur at double the Keplerian orbital velocity. Similarly, if there is variation 
in the orbital velocity of the prograde dust, some impacts into Europa's and 
Ganymede's surfaces could happen at significantly lower velocities. 



Fig. m 
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3. Structure of the dust model and calculation of input parameters. 

The dust cloud m odel used in these s tudies is based upon prev ious dust 



models developed by 



Kriiger et al. 



(120031 ) and 



Krivov et al 



(120031 ). We im- 



plement experimental and modelling impact data to reduce the number of 
unknown variables such as the slope of the ejecta size distribution and the 
size-velocity relation in ejected fragments. The relation ship between the in- 



flux a nd outflux is determined through the mass yield (IKoschny and Griinl . 



20011 ). The outflux populates the space around Europa or Ganymede, where 
part of the ejecta escapes (unbound dust), and part falls back to the surface 
(bound dust). In order to compare the contributions from different influx 
populations to the dust around Europa or Ganymede, outfluxes caused by 
each influx population were calculated separately and then combined. The 
model predicts the dust densities of the ejecta outflux at any altitude within 
Europa's or Ganymede's Hill sphere where these dust densities can be pre- 
sented as a size distributions of dust. Essentially, the model predicts the 
radial dust arising from the surfaces and can predict the size distribution of 
dust fragments in any 'onion-shell' around Europa or Ganymede. 
For each micrometeoroid impact, total ejecta mass is calculated (subsection 
3.1); then the total ejecta is transformed into ejecta size distribution (subsec- 
tion 3.2); and each ejecta fragment is allocated its corresponding ejecta speed, 
based on its size (subsection 3.3). Altogether combined, the total ejecta dust 
outflux from the surfaces of Europa and Ganymede are calculated (section 
4) and verified against existing data (section 5). 
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3.1. The ejecta mass yield and the total ejecta mass. 

The ejecta yield, Y, defined as the ratio of the total ejected mass, m c , 



)q. HI represents t 



l e effic iency 



200 If ). For 



to the mass of the impactor, m p , shown in 
of the material ejection in an impact event ( IKoschny and Grim , 
the same impact conditions, the excavated crater volume in ice-silica t e tar- 
get decreases with increasing silicate content , S (ILange and Ahrens 



Koschny and Grim , 



2001 



Hiraoka et al. 



19871 : 



20081 ). High S values correspond to 



the non-ice material and S=0 to pure ice. 



Y = m c /m p = V lfi (V hl0 o/V 1 ^ wo p t 2- b m b - 1 v. 



2b 
P 



(4) 



Vi o=6. 69x10 m J and V^ioo=10 m J are crater volumes made by 
an impact KE of 1 J, for S=0% and 5=100% silicate content, respectively, 
and pt the target density in gcm~ 3 , m p and v v are the mass and velocity of 
the impactor in g and kms~ l , respectively, and 6=1.23. The density of the 
target for different values of S is calculated using the linear mixing model 



( IKoschny and Griin . 



20011 1 shown in Eq. |5] (pj ce =0.927 gem 3 and silicate, 



Psii =2.8 gem 3 ). 



I/Pt 



S)/Pice + S/psU. 



(5) 



Europa's and Ganymede's surfaces are mostly covered with water ice, but 
contain some surface contaminants that are unevenly distributed. Europa 
is a bright icy moon with a young-looking surfa ce and a high albed o. In 



previous studies of the dust cloud around Europa Kriiger et al 



(120031 ) used 
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S'=0, but here we consider the non-icey part to have a non-negligible ef- 
fect on the ejecta, hence in Eq. [5] for Europa S'=0.1. This should account 
for some amounts of non-ice material observ e d in spectra of Europa's sur- 
face, both from Eart h (e.g. 



( IMcCord et al 



Spencer et al. 



(J2006D) and in Galileo flybys 



19981 ). From global geological mapping of Ganymede, it is 



estima ted that about 2/3 of its surface is covered with a dark non-ice ma- 



terial ( Patterson et al. 



20071 ). henc e we assume that £ =0.3 for Ganymede. 



This value is also in agreement with 



Krivov et al. 



(120031 ) who used the same 



value for their dust cloud around Ganymede calculations, but based the ap- 
proximation on albedo measurements. Both on Europa and Ganymede the 
dark non-ice material is mostly found in topographical lows. The mass yield, 
Y, and total ejected mass, m c , for all impactor dust grains that bombard 



Kriiger et al. 



(12003[ ) 



Europa and Ganymede were calculated using Eq. HI 
used a constant value F=10 4 , whereas in our model we use variable mass 
yields depending on impactor properties. 

3.2. The cumulative ejecta size distribution of the ejected fragments. 

According to the fragmentation law, the number of ejected fragments 
during a cratering event can be characterized by a power law function of 
the fragment's size. Therefore, the cumulative ejecta size distribution is 
represented in the form N = Ad~ B , where iV is the number of fragments 
larger than the fragment size, d. From impact experiments onto Europa sur- 



face analogue materials, 



Miljkovic et al. 



(120111 ) suggested that impact ejecta 
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from both pure ice and hydrated minerals fragment in a similar way, and 
that the fragmentation pattern is independent of impact conditions (such 
as impact angle). Assuming that these dist ributions are independ ent of the 



20091 ). they can 



impactor size, velocity and incidence angle (jMiljkovic et al. 
be applied to other impact events, such as dust clouds around Europa and 
Ganymede. In these previous experiments, the slop e of the cumulative eject a 



size distribution was 5=1.5 (jMiljkovic et al. 



20111 ) for both ice and non-ice 



materials. Assuming that these materials are suitable surface analogues, the 
results can be applied over the whole surface of Europa and Ganymede. 
The distribution coefficient, A, depends on the total impact ejecta mass (Eq. 
E]), but also on the target material properties: the largest ejected fragment, 
m b , that represents 1% of the total eje cta, M T , in each impact an d the ejecta 



fragments' size distribution slope, B ( lO'Keefe and Ahrens 



19851 ). 



A = M T 



1 - B 



Bm 



l-B 



(6) 



The cut-off for the smallest ejecta fragment size is set by the current labo- 
ratory detection capability (10~ 15 g, S. Barber, S. Sheridan, priv. comm.). 
The Io population was subsequently neglected as the largest ejecta fragment 
would be smaller than 10~ 15 g. 

3.3. The velocity distribution of ejected fragments. 

Each of the ejected fragments was assigned a velocity c alculated usin g 



Eq. [7] that corresponds to the ejecta size- velocity relation in 



Meloshl (119841 ) . 
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where v is the initial velocity of an ejected fragment, a t and p t are the surface 
tensile strength and density, respectively. 



v l V a t/pt = ci(d/d proj ) 



C2 



(7) 



The coefficient, ci, and exponent, C2, were derived from impact mo delling and 



impa ct experiments into water ice and non-ice brittle materials ( iMiljkovid . 



2010) and are equal to ci =0.06 and 



simulations are shown in 



c?=-1 .16. Details of numerical impact 



Miljkovid ( I2010I ). In summary, we used ANSYS 



AUTODYN-2D SPH (Smooth Particle Hydrodynamic) shock physics code 
to model impacts into ice and non-ice brittle rocky material. The ejecta 
fragments were resolved as clusters of SPH cells, from which their size dis- 
tribution relative to the size of the largest fragment was measured and their 
velocity recorded. For accuracy, the fragments were distributed in mass bins, 
each represented by X in Fig. [2j This size-velocity distribution w as verified 



agains t the experimental impact ejecta fragmentation study by iMiljkovic 



( 1201 ll ) to validate the numerical simulation results. By fitting Eq. [7] to the 



data in Fig. [2] we derived the coefficients c\ and c-i. Fig. [2] also shows the 



ejecta fragment size- velo c ity di st ribution from sim i 



Nakamura and Fuiiwara 


(1980 


); 


Nakamura et al. 



1991): 



(1994) 



Nakamuralfll993h : 



ar impact experiments by 



Fujiwara and Tsukamoto 



Fig. m 

From the perspective of impacting micrometeoroids, the top mm surface 



16 



layer at Europa and Gany mede is assumed to resemble snow ([Hansen and McCord 



2004 



SchenketaL 



20071 ) . so its strength can be expected to be lower than 
in the case of non-porous ice. We use a t =0 .3 MPa, corresponding to the 



approximate strength of snow (IPetrovic 



20031 ). The density, calculated from 



the linear mi xing model (Eg. El) was 0.99 gem 3 for Europa (10% silicate 



and 90% ice (IKrivov et al. 



2003)) and 1.82 gem' 3 for Ganymede (70% silica 



and 30% ice ( jMurchie et al. 



1992J) was 1.55); d and d pro j are the size of the 
ejected fragments and the size of the projectile, respectively. The combina- 
tion of ejecta fragment mass and corresponding ejecta velocity provides a 
good estimate of the dust outflux around Europa and Ganymede based on 
the modelled influx. 



4. Results. 

4-1. The ejecta outflux from Europa's and Ganymede's surface. 

For each impacting dust grain, the total ejecta mass, m c (m c ~ d~ 3 ), is 
calculated using Eq. HJ which is then transformed into an ejecta size dis- 
tribution in the form of N = Ad~ B . Coefficient A is calculated using Eq. 
[6] and B is the experimentally derived distribution slope. For each ejected 
fragment, the corresponding ejecta size-velocity relation (Eq. [7j) is applied 
and the ejected dust flux leaving the surface can be calculated. The sum of 
all ejected dust fragments created by the impactor influx, N, over the whole 
surface of a moon, corresponds to the total outflux, F out (Q), from the surface 
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(h=0) every second. 

Within the Hill sphere, the gravitational influence of the moon is the dom- 
inant force on ejected dust fragments. Using the energy conservation law 
(Eq. [8]) between the kinetic and gravitational potential energy of the ejected 
dust fragments, the speed of the ejecta fragments can be calculated at any 
altitude. 

rRm+h rh 

- / GMSdr/r = / d 3 v 2 dust (h)/2 (8) 

d and Vd us t{h) are the size and velocity of an ejected fragment at altitude h; 
G and R m are the gravitational constant and radius of Europa or Ganymede. 
By knowing the ejecta velocities at different altitudes, Vd us t(h) and the dust 
outflux ejected from the surface, F out (0), the spatial density at any altitude 
above the surface, N(h) (in m -3 ), can be calculated using Eq. [91 The spatial 
density was calculated using the ejected flux expanding through a sphere of 
a radius equal to the radial distance from the centre of the moon (R m + h). 

N(h) = F out (0)R 2 J((R m + hfv dust {h)) (9) 

Our calculations show of all five populations considered as influxes, only the 
asteroidal and halo populations are able to create micron size ejecta with 
sufficient speed to populate orbital altitudes around Europa and Ganymede. 
The other three (ISD, Io stream dust and ring dust) are found to create a 
small amount of ejecta that is either too small or too slow, so they have been 
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excluded from models. 

The cumulative spatial densities of ejected dust in the form of a mass dis- 
tribution at different altitudes above the surface were calculated using Eq. 
[9j These are shown in Fig. [3] for Europa and Fig. H] for Ganymede. The 
spatial density in both cases drops with radial distance. Ganymede's gravity 
is nearly twice the gravity of Europa which is reflected in the dust densities. 
We expect bound ejecta to eventually fall back onto the surface, effectively 
contributing to the dust population both on its upward and downward tra- 
jectory. We do not consider secondary ejecta due to the fact that bound 
ejecta moves at up to 2-3 kms~ l (which corresponds to the escape velocity). 
Figure 
Figure [4] 

If the kinetic energy {KE) of an ejected fragment is greater than the 
gravitational potential energy (PE), then the ejected fragment will escape 
the gravity of the moon becoming an unbound fragment. Otherwise, the 
bound fragments would fall back to the surface. The semi-major axis of a 
fragment's bound orbit can be estimated using a = —GM m /(2E), where 
E = KE — PE. 
Figures [5^ and b 
Figures [6^ and b 

Figs. [5^ and[6^ show the cumulative spatial density of ejected fragments 
larger than 10~ 15 g (which is considered to be the threshold for mass detec- 
tion in the laboratory) and Figs. [5b and[6b show fragments larger than 10~ n 
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g (the mass detection threshold in the Galileo DDS instrument) for Europa 
and Ganymede, respectively. Due to bound dust's upward and downward 
trajectory on which dust changes direction and falls back towards the sur- 
face, there is a prominent increase in dust density at around 500 km (0.3 Re) 
at Europa and a smaller increase around Ganymede at about 1000 km (0.4 
Rg)- Assuming that the top surface layer strength is similar on E uropa and 



Gany mede (0.3 MPa, which corresponds to the strength of snow (jPetrovid . 



20031 )). due to the higher gravity on Ganymede the distribution between 
bound and unbound fragments is different. On Europa 99% of ejected frag- 
ments belong to the bound population and will fall back to the surface after 
reaching the maximum height, whereas on Ganymede 80% of the ejected 
fragments are bound. A spacecraft in orbit around Europa or Ganymede 
could have orbital altitudes down to 200 km, hence it is important to study 
the dust populations at these altitudes to predict the dust count and size dis- 
tributions that can be collected and chemically analysed by a dust detector 
on such an orbiter. 



5. Verification of the model. 

This model was built on previous models upgraded with available impact 
experimental and modelling data. The model can be verified against the 
small amount of observational data collected by Galileo DDS (Dust Detector 
System) as only a handful of dust grains was measured in its short duration 
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flybys. 

During 8 Europa flybys, Galileo DDS collected about 50 grains (jKriiger et al 



20001 ) 



20031 ) and during 4 Ganymede flybys about 40 grains (jKriiger et al. 
that were believed to come from the dust clouds around the moons (N b s _). 
The available data were transformed into the cumulative spatial distribution, 
N = Ne or Nq in units of m~ 3 , by dividing the data by the spacecraft 
speed, v sp =2-9 kms -1 at Europa and 2-12 /cms -1 at Ganymede, the spin- 
averaged sensor area during the flybys (A =0.061-0.0235 m 2 ), the average 
time each flyby lasted nam ely about t=2 h and the number of flybys, Nfi y i, ys 



(jKriiger et al 



2000 . 



2003J) (Eq. [10 



N 



N, 



obs. 



(10) 



Vsp.tN flybys 

Fig. [7] shows the cumulative spatial density of dust around Europa, at 
600 km (0.4 Re) and 12000 km (8 Re) altitude , calculated by our mode l, 



Kriiger et al. 



( 2000 



20031 ). 



compared with the Galileo DDS data reported in 
where flybys were no closer than 600 km. Fig. |H] shows the cumulative spatial 
density of dust around Ganymede, at 600 km (0.2 Rq) and 20000 km (8 Rq) 
altitudes c ompared with the altitude s over which Galileo collected the data 



reported in lKriiger et al. 



(200 



o 



20031 ). Galileo flybys ranged between 0.1 and 



10 Rg at Ganymede. Figs. [7]and|8]show our dust model results a re consistent 



with the Galileo data at the corresponding detection threshold ( IGriin et al 



19921 ). apart from the small inconsistency among smaller fragments, which 
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are the most difficult to detect. The vertical error bar in Figs. [7J and [S] comes 
from the variations in the range of flyby velocities, the sensor area and flyby 
durations. The X-axis shows dust mass in bins. 
Figure [7] 
Figure [8] 

Figs. IH] and [TD] show the spatial dust density at different mass thresholds 
decreasing with distance, and is also compared with the Galileo DDS data 
measured at different altitudes from the surfaces of Europa and Ganymede, 
respectively. The calculated radial spatial density for m ass threshold corre - 



1992)) 



sponding to Galileo DDS sensitivity threshold (10 -11 g flGriin et all 
fits the Galileo data collected in Europa and Ganymede flybys. We are confi- 
dent that the predicted dust densities at lower altitudes, for which there is no 
observational data, are valid and can provide a confident dust estimate for a 
dust detector onboard an orbiting spacecraft around Europa or Ganymede. 
Figure [9] 
Figure [10] 

6. Discussion. Significance of input parameters accuracy. 

The most probable input values were used in the dust models presented 
here. However, the space dust environment is complex and possible input 
parameters could take a range of values. We investigated the effects a range 
of input parameter values would have on the outcome of the model. The 
variations of the ejecta mass yield and surface strength have the largest effect, 
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variation in the slopes of the cumulative ejecta fragment distribution or a 
change in the largest ejecta fragment affect the dust density to a lesser extent, 
whereas changing the rest of the input parameters (size- velocity distribution, 
micrometeoroid impact velocity, fraction of the non-ice content on the surface 
and micrometeoroid density) affect the model outcome the least. 
In our models the largest ejecta fragment was taken to be 1% of the total 
ejecta mass, M t , in an impact event. If that value was 0.01% M t or 10% M t , 
the dust density would increase or decrease, respectively, by a factor of 5. If 
the ejecta fragmentation slope B changes for ±20%, the dust density would 
also change by a factor of 5, moving the distribution more towards larger or 
smaller fragments, depending upon the slope value. 

The input parameter that has the greatest effect on the model outcome is 
the strength of the surface material, which is mainly ice. The strength of ice 
can vary significantly (it depends on temperature, grain size and porosity) 
from less than 0.3 MPa (corresponding to snow, which is the value us ed in 



our model), to 3 MPa for non-porous crystalline ice (jPetrovic 



20031 ) and 



up to 10 MPa if the ice contains a fraction of non-ice, silica-like, material 



(IHiraoka et al. 



20081 ). If the strength is varied between 0.01 and 10 MPa, 
this would create an order of magnitude (10 times) lower or higher dust 
density in the dust cloud, respectively. The strength of the surface material 
on Europa and Ganymede could also vary locally, which is another reason 
why it is important to emphasize the strength as the most dominant input 
parameter. The actual density measurements and chemical analysis of dust 
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in the clouds around Europa and Ganymede by a dust detector in orbit would 
also help narrow down the surface material strength. 

7. Conclusions. 

In this paper, we present an update model for determining the size and 
spatial density of dust around Europa and Ganymede at any altitude within 
the moons' Hill sphere, with particular emphasis on the role of micromete- 
oroid impacts ejecting surface material. The model computes influx distribu- 
tions from different sources and addresses each ejected fragment individually 
according to its size and corresponding speed, allowing the complete profile 
of the dust around Europa and Ganymede to be defined through calculation 
of the size, velocity and spatial distribution of the complete outflux. 
An important feature of this dust model is the ability to predict the dust 
around Europa for distances closer to the surface than any of the spacecraft 
flybys made so far as well as for the dust populations below the Galileo DDS 
detection threshold. We have verified these results by comparing our results 
with data collected by the Galileo DDS. It is found that the dust spatial den- 
sities close to the surface of Europa and Ganymede should be much higher 
than a simple extrapolation from the actual Galileo data. There should also 
be much more dust at sizes below the Galileo DDS detection threshold. 
This model can be applied to any other atmosphereless body in the Solar 
System, but above all, it is of great importance for future orbiter missions 
around Europa and Ganymede, such as recently selected JUICE mission, 
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ESA's next large space mission, to Jupiter System. This study provides sci- 
entific support for a dust detector /analyser payload for JUICE or any other 
future space mission to Europa and/or Ganymede. It also complements sur- 
face studies traditionally performed using remote sensing instruments. 
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Figure 1: Total impact flux (influx) onto the surfaces of Galilean satellites showing all 
dust populations considered as a source of dust around Europa and Ganymede. The in- 
terplanetary asteroidal (IDP-ast), ha lo (IDP - halo) and interstellar dust ( ISP) populations 
are presented as mass distributions ( Divinel . Il993t iLandgraf et all l2000h . Due to a lack 
of data on the Jupiter system dust (JSD), Io (JSD-Io) and ring (JSD-ring ) dust are pre- 



sente d as single points with estimated uncertainties (jKriiger et all l2004t iKrivov et al 
l2002h . The total influx includes the gravitational and magnetic field focusing and de- 
flection of fragments smaller than 10~ 15 g by the magnetic field from outside the Jovian 
system ( Colwell and Horanvi , 19961) . 
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Figure 2: Ejecta fragment velocity distribution normalized to \fotfpu where at is the 
impacted surface strength and pt surface density, calculated in ANSYS AUTODYN-2D 
impact simulations and comp ared to similar ejecta fr agmentation experiments by: [1 ] 
INakamura and Fuiiwaral (Il99ll) . [2] INakamura! (|l993l ). [3] lFuiiwara and Tsukamotol (|l980h . 
[4] iNakamura etlih ( 19941 ). pro vides a good fi t for C\ and C2 distribution coefficients, for 
the ejecta size- velocity relation (jMeloshl . Il984l ). 
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Figure 3: Cumulative spatial mass distribution of the ejected surface fragments at differ- 
ent altitudes above Europa's surface shows that the spatial dust density decreases with 
altitude and there are more smaller than larger fragments, which is both according to 
the fragmentation law and obeying the gravitational influence on ejecta fragments. Dust 
masses between 10 _15 g and 10 _11 g correspond to dust sized from 1 fim to 10 /im. 
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Figure 4: Cumulative spatial mass distribution of the ejected surface fragments at different 
altitudes above Ganymede's surface shows similar trends to Fig [3J Dust masses between 
10~ 15 g and 10~ n g correspond to dust sized from 1 /irn to 10 /im. 
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Figure 5: a and b. Cumulative spatial density of dust around Europa shown as a function 
of distance, r, from the surface (Re is radius of Europa), where dust masses are larger 
than 10~ 15 g and 10™ 11 g, respectively. Triangles show bound dust, squares show the 
escaping (unbound) dust and the solid line shows the total ejected dust fragments. Radial 
distance, r, is measured from the center of Europa in Re- Prominent increase in dust 
density at 0.3 Re is caused by bound dust reaching its maximum height, slowing down 
and reversing trajectory in that altitude range. 
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Figure 6: a and b. Cumulative spatial density of dust around Ganymede shown as a 
function of distance, r, from the surface (Rg is radius of Ganymede), where dust masses 
are larger than 10 -15 g and 10 -11 g, respectively. Triangles show bound dust, squares show 
the escaping (unbound) dust and the solid line shows the total ejected dust fragments. 
Radial distance, r, is measured from the center of Ganymede in Rq. An increase in dust 
density up to 0.4 Re is caused by bound dust reaching its maximum height, slowing down 
and reversing trajectory in that altitude range. 
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Figure 7: Cumulative spatial density of dust around Europa, at 600 k m (0.4 Re) and 12000 



km (8 Re) altitudes compared with the Galileo DDS data reported in ( Kriiger et al. . 2000, 



20031 ). where flybys were no closer than 600 km. 
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Figure 8: Cumulative spatial density of dust around Ganymede, at 600 km (0.2 Rq) and 



20000 km ( 8 Rg) altitudes compared with the Galileo DDS data reported in ([Kriiger et al 
2000L l2003h . where the flybys ranged between 0.1 and 10 Rg- 
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Figure 9: Radial spatial density of ejecta from Europa (measured from the centre of the 
moon), for dust masses larger than three different mass thresholds (10 -15 g, 10 -13 g and 
10~ n g), compared with G alileo data. Galileo DDS had a detection mass threshold at 
10" " 



g for this set of data (jGriin et all 119921 ) . which matches our model predictions. 
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Figure 10: Radial spatial density of ejecta from Ganymede (measured from the surface 
of the moon), for dust masses larger than three different thresholds (10 -15 g, 10 -13 g 
and 10 -11 g), compared with Galileo data. Galil eo DDS had a detection mass threshold 
at 10 -11 g for this set of data (jGrun et al which matches our model predictions, 

similarly to Fig. |9] 
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